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ABSTRACT 

It has been reported by different researchers that shear 
strain acceleration governs the machining parameters like, 
microhardness of chips, shear angle, tool wear, tool-chip 
interface temperature etc. Experimental and analytical studies 
have been made about dependence of tool-chip interface temper- 
ature on machining parameters during longitudinal turning. 
Experimental findings for above have also been reported during 
accelerated cutting. However, no literature is available to 
predict temperature during accelerated cutting (viz. facing and 
taper turning) analytically. 

In present work a general purpose finite element method 
(FEM) programme has been developed to calculate the average 
tool-chip interface temperature and temperature distribution 
in the tool-chip-work system during longitudinal turning as 
well as accelerated cutting .Suitable models for velocity field 
and heat generation in different zones have been used. 

The analytical (FEM) results have been compared acpinst 
experimental and analytical results of other researchers. A 
comparison of tool-chip interface temperature in accelerated 
cutting has been made with that in longitudinal turning to 
investigate the effect of shear strain acceleration on temper- 
ature attained during accelerated machining. 



It has been concluded that the average tool-chip 
interface temperature is highest in case of longitudinal 
turning and lowest in the case of facing. It is also 
observed that the temperature distribution pattern in all 
the three cases is approximately the same, except their 
magnitude. 



CHAPTER 1 


INTRODUCTION AND LITERATURE SURVEY 

1 .1 INTRODUCTION: 

One of the most important aspects of today’s metal 
cutting research is how to increase the productivity. The 
productivity itself is related to the life of cutting tool. 
Further, commonly used criteria for the evaluation of ease 
or difficulty with which any material can be machined 
(i.e. machinability) is tool life. Thus, to a certain extent, 
machinability evaluation helps in yielding higher productivity 
and developing new tool materials. Longitudinal turning is 
commonly used method for machinability evaluation. But, 
it suffers from inherent disadvantages of being slow, large 
amount of material consuming and expensive. Therefore, acceler- 
ated cutting tests viz. facing [1] and taper turning [2] were 
proposed for quick evaluation of the machinability. But, 

Jain and Pandey [3] reported that accelerated machinability 
test-results should be used with caution as they will be 
reliable only when shear strain acceleration is kept below a 
certain critical value, inspite of theoretical backing for 
these tests. 

The temperature attained during machining operation 
affects the properties of both, tool and work materials. In 
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other words, temperature is a critical factor which influences 
the shear angle, micro-hardness of chips, cutting forces, tool 
wear, surface finish etc. 

The following experimental techniques have been commonly 
used for evaluation of mean tool-chip interface temperature 
and temperature distributions in tool, chip and work surface [4] 0 

1 . Tool-work thermocouple, 

2. Radiation pyrometer, 

3. Colour paints. 

There are no simple and reliable experimental methods 
for determining the temperature field in the tool, chip and 
workpiece, even for the simple case of steady state orthogonal 
cutting. Therefore, analytical approaches although quite 
complicated, are relied on to make indepth analysis of thermal 
aspects of metal cutting process. 

Because of the high speed computation capabilities of 
moden age computers, numerical techniques like finite element 
methods (FEM) and finite difference methods (FDM) are being 
used to analyse the temperature field in machining. There are 
advantages of using FEM over FDM, such as a general purpose 
computer programme can be developed for a complex and changing 
geometry problems like the one encountered in metal cutting 
situation. By changing the input data only, it is possible to 
analyse the varieties of problems with different geometries 
and boundary conditions. It is not possible in case of FDM. 
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1 .2 REVIEW OF PREVIOUS WORK: 

Several investigators have studied, theoretically and 
experimentally, the temperature variation during orthogonal 
longitudinal turning. The average temperature rise in the 
cutting zone has been theoretically estimated by Zorev [5]* 
Reported estimation was based on the assumption that the total 
heat generated is distributed among the tool, chip and work- 
piece. Later on, theoretical analyses were also done by 
Weiner [6] and Rapier [7]. 

Finite element method (FEM) to determine the temperature 
around the cutting edge was first applied by Tay et al [ 8] . 

The concept of hyperbolic streamlines for obtaining the velocity 
fields and the strain rates in the primary shear deformation 
zone was used. Semi-empirical models have been used for calcu- 
lating the heat generation in the primary and secondary, zone. 

The FE solution was based on the variational formulation and 
quadrilateral elements were used for discretizing the solution 
domain consisting of a part of tool, chip and work material. 

Murarka et al [9 ] obtained the average chip-tool inter- 
face temperature using the Galerkin's approach of the finite 
element method for a wide range of cutting conditions. In 
his model, the shear strain-rate and velocity distributions 
in the primary and secondary deformation zones were determined 
experimentally using a grid technique. Empirical expressions 
for the flow stress, in the primary zone, in terms of the 
strain, strain-rate and temperature were used. Heat generated 
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at the interface of flank face and machined work surface was 
also accounted. 

Modifications in the model of Tay et al [8] were made 
by Stevenson et al [10]. The FE programme developed in the 
earlier work for calculating the temperature distributions in 
the chip and tool were extended in its range of application. 
Specifically, the programme no longer needed a flow field as 
input and it could accommodate a wide range of shear angle and 
contact lengths o The finite element temperature field results 
were compared with temperatures obtained using metallographic 
method. 

Balaji [ll] measured average tool-chip interface temper- 
ature while machining mild steel work using, both coated and 
uncoated WC tools. Using FEM, he computed average tool-chip 
interface temperature and temperature field in tool, chip and 
work material. He compared analytical results with his own 
experimental results. Simple linear model for velocity variation 
in the primary and secondary zone was assumed. 

Velocity components (Fig. 1 .1 ) u (along x direction) 
and v (along y-direction) were calculated as 

In the zone I Q , 

2V + V Sin a 
u = ^ 

V c Cos a 


v 


3 



5 



Fig. 1.1 Assumed deformation zone shape. 


v 
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In zone I Kf 


V + 2V C Sin a 


u 


V = 


2V C Cos a 


In zone II . 

O 


U 


^ 2V c + V int> 
< V c + V int> 


Sin a 
3 

Cos a 
3 


and In zone II b , 


u 


( V c + 2V int J 


Since 

3 


v = (V c + 2V. nt ) 


Cos ot 
3 


where. 


V int 


JL y ^ 1 + 8 “ 

3 c ° H 


The three noded triangular elements were used as discri 
tizing elements* 

The first attempt for measuring temperature during 
accelerated machining was made by Jain [12] . He measured 
tool-chip interface temperature (mV recording of thermo e.m.f.) 
during facing test, using H.S.S. tool and mild-steel as a work 
material. Jain and Pandey [13] determined temperature experi- 
mentally for facing as well as taper turning, using cast iron 
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as work material and H.S .S. as tool material. It was conclu- 
ded that strain acceleration plays an important role and it 
governs different parameters during accelerated cutting. 
Instantaneous temperature values obtained during facing and 
longitudinal turning for the same cutting conditions were 
found to be different. They concluded that accelerated tests 
for machinability evaluation are not reliable methods unless 
they are performed within the prescribed shear strain acceler- 
ation conditions, 

1 .3 OBJECTIVES OF PRESENT WORK: 

From the above literature survey, it is evident that 
no analytical work has been reported for temperature evaluation 
during accelerated cutting. Keeping this in view, finite 
element technique has been used to estimate the temperature 
distributions in the tool, chip and workpiece, during acceler- 
ated as well as longitudinal cutting situations. Temperature 
is also presented in the form of average tool-chip interface 
temperature for a range of cutting conditions. 

In the present model, suitable approximations for 
velocity field and heat generations in different zones have 
been proposed. By making an energy balance and assuming two 
dimensional heat transfer situation, the problem is analyzed 
using a general purpose finite element programme. 

The validity of the analytical model has been tested 
by comparing the analytical results with those obtained 
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experimentally by Jain [12] during facing and longitudinal 
turning. Finally, the effects of shear strain acceleration 
on temperature are investigated by comparing accelerated 
cutting results with those of conventional turning. 



CHAPTER 2 


THEORETICAL ANALYSIS 


Elaborate studies have been made in the field of 
mechanics of orthogonal metal cutting. Here, theoretical 
analysis for the case of accelerated cutting, viz. taper 
turning and facing operations, in which cutting speed varies 
continuously during the cut, has also been included. 


2.1 SHEAR ANGLE, SHEAR AND FRICTIONAL FORCE AND SHEAR FLOW 
STRESS: 

In an attempt to study stress-strain properties 
during machining process, Kececioglu [14] obtained photo- 
micrographs of the plastic zone, using quick stop device 
(QSD). Based on these photomicrographs, he assumed that 
plastic zone could be represented by a parallel sided shear 
zone, similar to that shown in Fig. (2.1). 

The line AB is the centre line of actual shear 
zone DCEF and known as shear plane (Fig. 2.1). The position 
of shear plane is defined by 0, known as the shear angle 
and it can be evaluated (Fig. 2.3), from the equation (2.1) as 


tan 0 


r Cos cc 
1 - r Since 


(5>.1) 


where , 


r is chip thickness ratio given by. 
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AREA d'c'eV ezza Plastic deformation zone obtained from photomicrographs. 
AREA DCE.F 6sa Approximation of o'c'e'f' in to a parallel sided zone. 

Fig. 2.1 Chip formation and deformation during orthogonal cutting 




Fig 2.3 
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r = 




( 2 . 2 ) 


where. 


Lj and lg are lengths of undeformed and deformed chip 
respectively 

f = feed rate (undeformed chip thickness) 
t = deformed chip thickness. 

Forces can be evaluated using Ernest and Merchant's model of 
thin shear zone [15] . Merchant's force circle diagram is 
shown in the Fig. (2.2). 

Mean shear stress, r , is given by 

( F Cos 0 - F. Sin 0) Sin 0 

■r = — r-r (2.3) 


where, 

d = depth of cut. 

2.2 SHEAR STRAIN, STRAIN RATE: 

Shear strain is defined as the deformation per unit 
length, i.e. if a rectangular block of height 'a' slips by 
an amount 'b' (Fig. 2.4) then, strain. 


v 


b 

a* 


(2.4) 




Fig. 2 A Concept of shear strain 



Fig. 2.5 Elements of chip in strained state. 



A 


Fig. 2.6 Enlarged view of area ABC of fig. 2.5 
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The chips can be considered as composed of a large number 
of elements of thickness 'a 1 and then each element suffers 
a displacement ' b’ after crossing the shear zone* Fig, (2. 5). 
Then, it is clear from Fig, (2,6) that 


Cos f* 

Sin 0 o Cos( 0 “-a) 


(2.5) 


If 1 dt’ is the time taken by the metal to move through 
a distance in PSDZ and V is cutting speed, then strain rate 


v 


dv 

cFE 


( 2 , 6 ) 


V Cos oc 

2^ Cos(0 - a ) 


(2.7) 


where. 


•d.a i = thickness of PSDZ 

4 


Hence, 

♦ 

V 


V Cos a 
D'S .Co“sT0- a)" 


DS . 


( 2 . 8 ) 


Jain and Pandey [16] have shown analytically that 
the mean width of PSDZ can be given by 


V 2 f Sin( 90 + a - 0) 

4 Sin 0.Sin(45-a+ 0) 


(2.9) 
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2.3 SHEAR STRAIN ACCELERATION DURING TAPER TURNING AND 
FACING OPERATIONS: 

During taper turning, cutting speed increases, as 
cut proceeds from smaller end towards larger end, as shown 
in Fig. (2.7). The speed at any instant ' t' is given by 

v t “ Tffik [ d b + 2 tan e t .Nf-t] (7.10) 

For facing, the cutting speed at any time * t' after 
the cutting starts at inner dia D^ Fig. (2,8), can be written 
as 


V 


f 


TCN 

TooS 



+ 2Nf t 3 


( 2 . 11 ) 


If the variations in *0* for taper turning and facing can be 
represented in the following forms, assuming constant depth 
of cut, then 


0 * - a, n" 11 . e " 12 . f" 13 . v^ 14 


m 


11 


m 1 3 ,, m 14 


0 f = B 1 N 1 1 . f , ' J . V f 


( 2 . 12 ) 

(2.13) 


where the constants A 1 , n^, n 12 , n 13 , n 14 , B 1 , ra^, m 13 , m l4 
can be found experimentally. A typical set of values of these 
constants are given in Table (2.1). 
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TABLE 2.1 : Constants Used For Calculation of 0 in the 

n-j -j n 0 ni o 
expression of form, 0 = A N O^ 1 ^ f J 

V* 1 4 (. -AWtd ) 


Taper Turning 

Facing 

A = 30.5395 

B = 82.2896 

n^= 0.0141 

m 1 -0.1729 

n 1 2 = 0*056 

m 12 = 0.0000 

n 13 = 0.2352 

3 = 0 . 0626 

n 14 = -0.1418 

m 14 = 0.1404 


When cut proceeds, all quantities in R.H.S. of 
eqns. (2.12 and 2.13) except 'V' remain same, so one can 
write the above equations in the form 


0 



(2.14) 


where , 


'0 


is a constant given by 


a 0 = vn" 11 . e t " 12 . f" 13 


(for taper turning) 


and. 


(2.15) 
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,n 


m 1 1 ® 1 3 


B 1 . N ' 1 . f 


n 


V“ = V, 


14 


V' 


n 


m 


14 


(for facing) 

(for taper turning) 
(for facing) 


Strain rate from eqns . (2.8) and (2.9) becomes 


V 


4V Cosof. Sin 0 . Sin(45 + 0 -a) 
42 Cos 2 (0 -a) 


Strain acceleration is given by. 


d£ 

dt 


Differentiating (2.19) logrithmically vdth respect to ' t* 
we have 



where. 


V = 


dV 




0 


dg 

dx 


(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


(2.23) 
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Differentiating (2,14) w.r.t. * t’, we have 

b = T55 A 0 • n Vn ”’ V (2.24) 

n . V , v 

= n $ • V (2„25) 

Since 0 is in degree in eqn, (2,14), 

From ecps ,( 2,21 - 2,25), 

v - v ^ [l + :j-g^ 0 n {Cot 0 + Cot (45+0-a)+2 tan(0-a)}] 

(2,26) 

Hence, strain acceleration for taper turning can be evaluated 
as, 

^ *t ^*t v" ^ ^ ^ TOO ^*t n l 4 ^ 0^ Cot( 45 + 0^ *■*&) 

"t 

+ 2 tan( 0 t - a) }] (2,27) 

where , 

can be found by differentiating ( 2 , 10 ) w.r.t. 

1 1 ' and given as 

2jiN 2 f 

Vt = Yoqo tan » t (2,28) 

and 

V t , n 14 , 0 t , y t are given by ec[ns,(2.10, 2.12, 2.19), 
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Similarly, strain acceleration for facing is given by, 

* • v f V- % 

v f = [1 + ^ 0 f . m 14 C Cot 0 f + Cot( 45+0 f - a ) 

+ 2 tan (0 f -a) 1] ( 

where, 

• 2 nN 2 f , 

v f = 1000 ^ 

(obtained from eqn c 2.1l) 

• © 

V f ,0 f , m 14 , v f are given by eqns. (2,11, 2.13, 2.19). 

In the work reported in this thesis, eqns. (2.10-2.13 
2.19, 2.27-2.29). have been used for computation purposes. 


2.29) 


2.30) 



CHAPTER 3 


CO-ORDINATE, VELOCITY FIELD AND HEAT GENERATION 

3.1 DIMENSIONLESS CO-ORDINATES: 

Unlike in longitudinal turning, cutting speed changes 
continuously in accelerated cutting as the cutting progresses 
with time. This change in cutting speed alters the values of 
different cutting parameters like shear angle, thickness of 
PSDZ, chip contact length, forces etc. [17,16,8] . Thus, 
in turn, the configuration of the solution domain changes 
with time (cutting speed), and co-ordinates of any point in 
the domain are functions of time. To take into account the 
grid size variation, the domain is divided into several zones 
Fig. (3.2), and the co-ordinates of a point in any of these 
zones are non-dimensionalized with respect to the two adjacent 
sides of a parallelogram, either made by the sides of that 
particular zone or by enlarged lines containing that zone 
sides. The co-ordinates in the non-dimensional form are 
kept the same with time, while the actual zone configuration 
and the orientation may change. The current configuration 
at any time ' t* can be mapped into the initial configuration 
(at time t=0), by suitably transforming it through translation, 
rotation and scaling. 
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Thus, at any time a transformation matrix can be 
found which when multiplied with non-dimensional co-ordi- 
nates gives the co-ordinates of any point at that instant. 

Mathematically speaking, 

CX p (t)} = [ TR( t) ] {X pd } (3.1) 

where, 

CXp(t)! = co-ordinates of a point at time ’ t' 

[TR(t)] = co-ordinate transformation matrix at 

time «t' 

{Xpd? = non-dimensional co-ordinates of point, 
given at time t = 0 (remains fixed) . 

Such a transformation for a parallelogram-shaped 
zone is shown in Fig. (3.1). 

In the figure, ABCD is the zone configuration at any 
time ' t’. It can be described by dimensions 1^, 1^ and the 
angles, 0^, © 2 * and y A are the co-ordinates of the corner 
point, A. 

Let (xp, Yp) be the physical co-ordinates of a point 
P, within the zone ABCD. In the transformed space, x pd and 
yp d are the dimensionless co-ordinates. 

Then, 



22 



Fig. 3.1 Parallelogram shaped zone ABCD. 



~ cosec 0 


Fig. 3.2 Division of solution domain (twelve zones) 
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X = £9 

Pd 1 1 

(3.2) 

QP 

Y *> - % 

(3.3) 

OP = (AQ + QP) + OA 

(3.4) 


From eqns.(3.2-3.4) and Fig. (3.1), 

xp 1 + y p J = (x pd .l 1 Cos © 1 + y pd 1 2 Cos© 2 + A x )l 


+ Cxpd* 1 -] Sin© 1 + y pd 1 2 « Sin© 2 + A y )J (3.4a) 


Above equation can be written in the matrix form as 


x p 

YP 

1 


} 

Cos0«| 


A x' 

r x Pdi 

L " 



\ 

) 

f 

fej 1 - -C03>©^ 

r pd 

J 

0 

0 

1 



(3.5) 


Figure (3.2) shows the geometry parameters of different zones. 
The parameters have also been tabulated in Table 3.1 . These 
parameters have been used for calculating transformation matrix. 
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Continued 


TABLE :3.1 ( Continued’ 



Co-ordinates (dimensionless) for a point in zone II is given by non-dimensionalizing 
sides hd and he s 
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3.2 ESTIMATION OF MACHINING ZONE DIMENSIONS: 

The length (1) of the primary shear deformation zone 
(PSDZ) along the shear plane is given by 

1 = f/Sin0 (3.6) 

Secondary zone is assumed to be triangular in shape. The 
maximum thickness of secondary shear zone (6t 2 ) has been 
reported by Tay et al [8] as: 

6t 2 = 0.05 t c (3.7) 

In the present analysis, the secondary plastic zone 
is assumed to extend over the entire contact length, neglecting 
the elastic contact region which occurs just before the chip 
leaves the rake face. A good estimate of the chip contact 
length can be obtained by assuming that the point G at which 
plastic contact ends is located such that AG is parallel to 
the resultant force, F R Fig. (3.3a) 0 

Then, from geometry, 

H _ 1 

3Tin§ Sin( 90 - (a - 0 + 0) ) 

from eqns. (3.6) and (3.8) 

f Sine 

Si n0 . Cos (a - $ + ©J 


H 


(3.9) 



Chip 



Fig. 3.3(a) Force diagram 


Fig. 3.3(b) Hyperbolic model of 
a streamline. 
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3*3 VELOCITY FIELD: 

3 .3. (a) For Uncut Metal and Chip: 



The 

velocity of uncut metal is the cutting speed (V) 

itself. 

So the 

velocity components for the uncut metal are 

• 

• 


u 

= 

V 

(3 


V 

= 

0 

(3 

where. 






u 

ss 

component of velocity along x-direction 



V 

= 

component of velocity along y-direction. 



The velocity components of fully formed chip are: 



u 

as 

V c .Sina 

(3. 


V 

5=5 

V c . Cos a 

(3. 


where, 

v c = V.r . 

V c is chip velocity and r is chip thickness ratio. 


3.3(b) Primary Shear Deformation Zone (PSDZ): 

The hyperbolic streamline Fig.(3.3b) concept developed 
by Tay et al [8] is used for calculating the primary zone 
velocity. Strain-rate and heat generation-rate. According to 
this model, a material particle is considered to flow along a 
hyperbolic streamline curve in the PSDZ. The model is based 
on the following assumptions: 
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(i) The direction of the maximum shear strain-rate is 
parallel to the shear plane. 

(ii) The overall change in velocity through PSDZ can be 
represented with sufficient accuracy by the velocity 
diagram Fig. (3.4). 

(iii) The distribution of maximum shear-strain rate along 
the normal to AB Fig. (3.3a) is symmetrical about 
AB with the maximum value occurring at AB. 

(iv) The strain-rate does not vary in a direction parallel 
to AB. Therefore, all streamlines passing through 
PSDZ can be taken to be similar in shape. 

All the above said assumptions had been claimed to be 
verified experimentally by the authors of the model. Envoking 
these approximations, they have derived the equation (3.14) of 
the hyperbolic streamline 

y (y tana - x) = a (3.14) 


where. 


a = constant. 


Using velocity diagram. Fig. (3.4), velocity Vp 
point on a streamline is given by 

V Sin 0 

V p * Sin(szS + 


at a 


(3.15) 


where P is the angle made by the streamline with the cutting 
direction. It is obtained by differentiating eq. (3.14) in 
the form. 
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tan (3 


dx 


2y tan a - x 


(3.1* 


Velocity components are 


u 


= Vp . Cos p 
v = V ft . Sinp 


(3.17 

(3.18 


The component of velocity parallel to AB, will be 


V {30 = Vj3 . Cos(0 +P) 


The strain rate is then calculated as: 


O 

V 


d V, 


M 


dz 


(3.19 


(3.20 


where z is the distance of the concerned point from the mean 
plane AB, Fig. (3.5), and it can be seen that 

z = x Sin0 + y Cos 0 . ( 3 .21] 

From eqns. (3.14-3.21), we get, 

• 2a V Sin0 

^ s o/o (3.22) 

[z + 4a Sin z 0(Tan a + Cotft)]-*' z 

The constant ’a 1 of the eqn. (3.14) can be etensimted 
as follows. 

The shear velocity reaches the value V within the chip 

s ~ 

from zero within the uncut metal. Thus, average work done per 
unit time in shearing the metal can be approximated as 
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W 

s 


F 


s 


V 

s 


2 


(3.23) 


as: 


The work done per unit volume, w , can now be written 

b 


W 


s 


F s 

dat Bs) 



<V> 


T 


V 

2Cd§7 


(3.24) 


Based on the assumption of a constant yield shear stress 
in the PSDZ, the average shear strain rate can be shown to be 


V 


av* 


w 


v s 

T • 2(DS) 


or 


w 

av # 


V s 

2(DS7 


(3.25) 


Referring Fig. (3.5) and eqn.(3.22), the average strain rate 
is also given by the expression: 


. 

v 


av 


DS 



o 


_2 

DS 


DS 

2 

I 

o 


dz 



J 

+ K 2 ) 3//2 


dz 


(3.26) 


where, 

J = 2a V Sin0 

and, 

K 2 * 4a Sin 2 0 (tana+ Cot0) 


(3.27) 


(3.28) 
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Now, 

J = V Cos a 
2 Cos( 0- a) 

V 

= (using velocity diagram 3.4) 


Integrating eqn.(3.26) with the substitution 
z = K tan©, we get 


v 

av 


2 _ 

DS 


DS 

2 j 

/ . Cos© d© 

z=0 K 


2_ 

DS # 



(f ) 2 + 



which, using eqn.(3*29) becomes 


V = 2 V If?T ' I 5 

Equating eqns.(3.25) and (3.30), 

K 2 = | (DS) 2 . 

So, strain rate can be written as 

^ ( z 2 +k 2 ) 3/2 

where, 

J, K and z are given by the eqns. (3.27), (3.28) 


(3.29) 


(3.30) 


( 3 .31 ) 


(3.32) 


and (3.21) respectively. 



33 


3, 3(c) In Secondary Plastic Zone: 

By analysing the deformation along the tool-chip 
interface obtained using printed grids with quick stop device 
( QSD) , Tay et al [8] found that the sliding velocity of 
the chip at the interface could be given by the equation 


V 


int 




(8.33) 


where , 

S is the distance of the point concerned with the cutting 
edge measured along the interface Fig. (1.1 ). 

Velocity in the secondary zone is linearly interpolated 
from at interface to V c at the fully developed chip. 

The velocity field thus obtained is used in turn for 
calculating the strain rate and the heat generation rate within 
the secondary zone. 

3.4 CALCULATION OF HEAT GENERATION: 

In the primary zone, heat generation occurs due to 
plastic work only. In the secondary zone it occurs both due to 
plastic work and interfacial friction. 

Assuming the shear stress to be uniform, the rate of 
• heat generation per unit volume in the primary zone is calculated 
as 
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(3.34) 


The heat generation rate per unit volume in secondary 
zone [8] due to plastic work can be assumed to vary linearly 
from QsZintp i n *t er f ace to zero within the fully formed chip* 
The interfacial value ^SZintp is given by 


^SZintp 


10 F f V 


3 St 2 (3H - b)d 


(3.35) 


where. 


6 t£. tan (0 -a) 


(3.36) 


In secondary zone, frictional heat generation rate per 
unit volume is obtained by distributing total frictional energy 
over a small strip of average width. Ay, and length, H, in the 


form: 


F f V c 
3Hd A y 


1 + 


(3.37) 



CHAPTER 4 


FINITE ELEMENT ANALYSIS 

INTRODUCTION [18 ] : 

Finite element method (FEM) is a powerful numerical 
technique for solving many problems of engineering importance, 
such as the heat conduction in solids, stress analysis for 
design etc. The FEM can be applied without much difficulty 
to problems with irregular geometries, complex boundary 
conditions and non-linear processes. The advantage of the FE 
procedure over other numerical techniques such as the finite 
difference method (FDM), is that general purpose computer 
programmes for an application may be developed through FEM, 
which can be applied to a wide variety of boundary conditions 
and geometries by just changing the input data. Due to these 
reasons, the FE solution technique was employed in the present 
work for predicting the temperature field during accelerated 
mac hining. 

4,1 DESCRIPTION OF FE TECHNIQUE: 

The solution steps in FE can be summarized as [18]: 

(a) obtain the equivalent integral equations from 

governing differential equations. The boundary 
conditions would be transformed into appropriate 
boundary# integrals. 



36 


(b) Divide the solution domain into subdivisions or 
elements* (Figure (4.1) shows such subdivisions 
to the problem analyzed in the present work). 

Place chosen number of nodes within each element 
and obtain approximating functions for describing 
the field variable, within each element, in terms 
of nodal values. 

(c) Find the elemental matrices for the nodal values 
within each element. 

(d) The elemental matrices are then assembledfor all 
the elements, giving rise to global matrix 
equations for all the nodal values in the entire 
solution domain. These global matrix equations 
are solved to get the value of field variable at 
all the nodes. 

Variational or Galerkin’s procedure can be followed 
for converting the governing partial differential equations 
into equivalent integral equations. 

The Galerkin's approach is used in the present work. 
Consider the differential equation of the form, 

L(0) = 0 (4.1) 

where 0 is the exact solution to the differential equation. 

* 

However, approximation solution, 0, will leave behind a 
residue . 
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That is, 

L(0*) = R ( 4*2) 

where R is the residue which may be a function of the spatial 
co-ordinates. In Galerkin’s procedure, the residue, weighted 
approximately by suitable weighting functions, is minimised 
over the whole solution domain D in an integral sense. In 
mathematical form, this becomes 

f D W.R dV = W j[ L(0*) dV = 0, for all i (4.3) 

where are the weighting functions. The approximation solution 
* 

0 may be written as a piecewise continuous profile for each 
element as: 


0*(®) = “ 0^ (4,. 

i=1 1 1 

where m is the number of nodes per element. In the above 

equation, 1 are the approximating functions for representing 

0 within the element (,e; and these are called the "shape 

functions" of the element (e). The nodal values of 0 for 

the element (e) are denoted by 0^ e ^ . In Galerkin’s procedure, 

the weighting functions, W. are taken to be the same as the 

( e) 

shape functions ' within each element. 

The governing differential equation for heat transfer 


during machining can be written in the form. 



38 


k + 3 %) 

ax ay 


pc p (u lx + v ^ + Q 


= p c 


3T 
p at 


( 4 .4a ) 


Using an approximate temperature T in eqn.(4.4a), we have. 


a 2 T a 2 ! a = a = 

L( T ) = k ( — 5 + — y) -P c (u |~ + v |“) 

ax ay p ax ay 


+ Q -PC P || = R 


(4.5) 


For the Galerkin’s formulation. 


t(e) 


9- 

a t 


d 2 ! 


2 //N^ '[k ( — 2 + 2) ** PC n 


e=1 


ax 


ay 


al 

ax 


+ v |I) 
ay 


a T 


+ Q “pc p ;-^]dx dy 


0 


(4.6) 


i = 1,m, where m is the number of nodes per element 
and E is the total number of internal elements within the domain 
D . Now, 

m 

T = 2 N^ e ^ T^. e ) for element (e) 

j=1 ^ J 

= lN (e) J t4 e) ) (4.7a) 

( Q \ ( 0 \ 

where LN X y J and {T^ '} are the vectors containing the shape 
functions and nodal temperature values for element (e). For 
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representing T over the whole solution domain. The above 
equation may be written as, 

T = LNJ{T n } (4.7b) 


whereLNJis the overall shape function vector containing the 
piecewise continuous approximating functions and £ T n ) is 
the total nodal temperature vector. Using eqn. (4,7a) in 
eqn, (4.6), we get 


E 

2 

e=1 


m 


//(e) £ kN i 


(e) 


3 2 n^ 

(— *• 


4 * 





(e) 3N[ e) 3N (e) 

~ K i P c p (u — 


JL + v ^J_) T (e) 


+ Q - 


,(«> ill 


(e) 


i - - N i c P tt 


+ T 


(e) 


3 N 


(e) 


j 3 t 


4—) ] dx dy 


0 


( 4 . 8a ) 


j = 1,m 


The integral in the above form would require the 
continuity of shape functions and its derivatives at all 
the interface regions, since second order derivatives of Nj 
occur in the equation (4.8a), It is useful to convert the 
above equation by making use of Green's theorem, so that 
lower order shape functions can be used for approximating the 
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temperature. The Green’s theorem states that for two scalar 
functions f.j and f 2 . 


/// v •( vfn) dV = 
V 1 z 


& n 
S “ 


* f 1 v f 2 ds 


or 


tn [ 


• vf2 + f-j 



v 


f, 

3^ 


] dV 
ds 


(4.9) 


(4.9a) 


using the above theorem for the functions N. and T. 

1 * 


/// k N, . v T dV 
D 1 


1 k N i If ds 

n 


/// U r 3N j 


V 


k C 


3T 

3x 0 ax 


3 A 31 

ay *3y 


a n 


a t 


Jt dx dydz 


(4.10) 


Combining the results of eqns.(4.10) and (4.8a) in 
2-D co-ordinates, we get 

3 N 

V LINIJ 

Tx 


u k r 3 N i 3 [n] 

D k L T* • "fx ay 


+ -i - 4NL ]iTn , d!t dy 


a y J n' 


+ k Nj 


StMl 


a x 


n 


dl + D 3 N, dx dy 


-•y pc p N i£ u + v tT n ) dx dy 
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" S D pc p N i[L N J*3r } + ^f^-CT n }]dxdy 

= 0 


On the boundary (boundary conditions used for present work 
is shown in Fig. (4.2) of the solution domain D, we get 


f k N, ~~ — dl 


i 3x 


n 


s N i t - k Ir 3 dl 

Oq 1 9X n 


S h C N i - k dl 


l N. q dl - / Ni h (T-T ) dl 

o q i o h o 


( N. q dl + l N. h T dl 
S q i ^ i_ o 


L N.h T dl 
S h i 


= - jC N, q dl + £ N.h T rt dl 

Oq i 1 ° 

- ^ Ni.hLNJ{T n } dl 


Therefore, 

g E N ± dA 



3LML + L!k HnJ ) 

3 x 3y 3 y ' 


+ P<=p 


^ (u 


§JiLL + v 2JLH1 

ax ay 


*1*1)] 


(4.11) 


(4.12) 
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r\ 



FIG. 40 FINITE ELEMENT MESH 
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Fig. 4»2 Problem region, showing the thermal 
boundary conditions. 



lb} b- s domain 



F, 9- 4.4 A typical boundary C Line) element 
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dx dy {T n ! + " p c p H ± {. N J dx dy l } 

S h \ L N J tt n ! dl ** s N i q dl 

Q 

+ S h N i h T o dl + D* Q % dx dy =0 (4,13) 

The eqn.(4,13) can be written down in the matrix form 
as, 

3T 

[H J CT n } + [Mj C-^} » {F} (4.14) 


where, for the matrices [H] and [M] , and vector ?F? , the 
elements are given by 

H if = 4^ + H 1 S) . (4.14a) 

i 3 

where. 


H\ e ^ 

" (e) hN i N J 

dl 


(4.15) 

r/( 6 ) 

K ij 

3N 

= (i) t k (?r 


3 N 

4- mmmmmm 

ay * 

3N 

mrnJm ) 

ay ' 


a n 1 

+pc P (u N i r * 

3 N . 

+ ^ dx dy 

+ V Nj. 

3 N. 

ay 

(4.15a) 


This can be written as 

[H] = E [H ( 

e=1 


( 4 ,15b) 
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s [K (e) ] + [H (e). 

e=1 1 


e=1 


II 


ij - (e) p c p N i N j dx d y 


(4.15c) 

(4.16a) 


This can be written as 


[M] 


Z [M (e) ] 



e=1 



F (e) _ 

II 

(e) 

Q Nj dx dy - 

( e) N i dl 


+ 

,£\ N. h T 

(e) i o 

dl 

which can be 

written as 


{FI = 

E 

2 

{q C e )j + 2 E q 

«4 #) » 


e=1 

e=1 



E h 

4 2 

e=1 

C4 e) ! 


(4.16b) 


(4.17a) 


(4.17b) 


where. 


i = 1,2, ,,..m 

j = 1,2, ,...m 


and E . E. are respectively the number of boundary elements 

4 h 

having applied heat flux and convective boundary conditions. 
In eqn. (4.17b) 


•(e) // • 

CQ 1 = ( e ) Q % dx dy 


(4.18a) 
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*H ( 2 e) ! 


N q dl 


(e) "i 


(e) N i h T o dl 


(4.18b) 

(4.18c) 


4.2 SHAPE FUNCTIONS AND THEIR LOCAL DERIVATIVES: 

In the present work, 8-noded quadrilateral isoparametric 
elements are used for spatial discretization of solution domain 
Fig. (4.3). 

Let (x^t Y^) be the cartesion co-ordinates of i node 
in x-y domain and after mapping the element in r-s domain 
(r^, s i^ are its natural co-ordinates. Then natural and 
cartesian co-ordinates are related by 



N — |jN.j ^ ••• ^8_j 


(4.1 8d) 


(4.1 8e) 


(4.1 8f) 


(4.1 8g) 


Nj’s are functions of r,s and given by (4.21a - 4.21 g) 
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0 = 1 0 0 0 .... 0 J 


( 4 . 1 8h ) 


Temperature within each element is in the form 

T = a Q + a 1 x + a 2 y + a 3 x 2 + a 4 y 2 + a & xy 

2 • 2 
+ a 6 x y + x ? x y c 


(4.19) 


In terms of the nodal values of temperature and the 
shape functions , we can write eqn. (4,19) as 


8 

T = £ N, T. 

i=1 1 X 

where, 

N 1 = 0.25 (-1 + rs + r 2 - r 2 s - rs 2 ) 

N 2 = 0.5 (1 - s - r 2 + r 2 s) 

= 0.25 (-1 - rs + r 2 + s 2 - r 2 s + rs 2 ) 

N 4 = 0.5 (1 + r - s 2 - rs 2 ) 

= 0.25 (—1 + rs + r 2 + s 2 + r 2 s + rs 2 ) 

N 6 = 0.5 (1 + s - r 2 - r 2 s) 

N 7 = 0.25 (-1 - rs + r 2 + r 2 s - rs 2 ) 

N 8 = 0.5 (1 - r - s 2 + rs 2 ) 


(4.20) 


(4.21a) 
(4.21b) 
(4.21c) 
( 4.21 d) 
( 4.21 e) 
( 4.21 f) 
(4.21 g) 
( 4.21h) 


Differentiating the expressions (4.21a - 4.21 h) partially 
with respect to r and s, we get local derivatives of shape 
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functions as. 


3N 1 

3 r " 



0,25 (s + 2r - 2rs - s 2 ) 

(4.22a) 

3N 2 

Sir" = 

-r + rs 

(4.22b) 

3 NL 



3 r 

0.25 (-s + 2r - 2rs + s 2 ) 

(4.22c) 

3^4 _ 

3 r = 

0.5 (1 - rs) 

( 4.22d) 

3N. 5 

51T = 

0.25 (s + 2r + 2rs + s 2 ) 

( 4,22e) 

aN 6 



o ^ 

ix 

-r - rs 

( 4.22f ) 

3N^ 

cflT 

0.25 (-s + 2r + 2rs - s 2 ) 

(4.22g) 

3Ng _ 

3 r " 

0.5 (-1 + s 2 ) 

( 4.22h) 

3R, 

3s 

0.25 (r + 2s - r 2 - 2rs) 

(4.2.3a) 

3N 2 

3s “ 

0.5 (-1 + r 2 ) 

(4.23b) 

3N 3 _ 

3 s 

0,25 (~r + 2s - r 2 + 2rs) 

(4.23c) 

3^4 = 
3s 

-s - rs 

(4.23d) 

^5 = 

3S " 

0.25 (r + 2s + r 2 + 2rs) 

( 4.23e) 

3 s " 

0.5 (1 - r 2 ) 

(4.23f) 

3N ? . 

3 "s *“ 

0.25 (-r + 2s + r 2 - 2rs) 

(4.23g) 

^8 
as * 

— s + rs 

(4.23h) 
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4.3 GLOBAL DERIVATIVES AND NUMERICAL INTEGRATION: 

(a) Global Derivatives: If G is a scalar function of 
the natural co-ordinates r and s, its derivatives with respect 
to x and y can be given by 


3G 

3G 

3x . 

3G 

iz 

3r 

(4.24a) 

3 r 

3"x ° 

S'? + 

ay e 

9G 

as 

= is 

” ax * 

is + 

3s 

3G 
ay • 

ix 

3s 

(4.24b) 
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functions can be written as 


1LN-L~ 


3-1. Mi 7 

3 X 


ar 


= [J]” 1 


31 N J 


3 1nJ 

U 3 y 


3s J 


(b) Numerical Integration: In order to get the elemental 

matrices, the area integrals in the eqns. (4*1 6a), (4.17a), 

(4.15) etc. have to be evaluated, which are in the form 
3 N* 3N , 

(l) f (N i* 3x > ay ••••)* dy • (*.71) 

The quadrilateral region in the physical co-ordinates 

(x,y) has been mapped into square region in the natural 

9 3 

co-ordinates, and the global derivatives g-j * — as a 

functions of natural co-ordinates are known, therefore, 

these integrals can be written in the form 

1 1 

I ~ J\ -i g( r »s) dr.ds (4.28) 


where, 

dr ds is given by 

dx dy = det. [J] dr.ds (4.28a) 

Now integrals in the above form can be evaluated 
numerically by using Gauss quadrature formula, taking -9 Gauss 
points per element 
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I = 

L £ 

(4.29) 

where. 



the number of Gauss points in one direction. 


n 

= 3 


and W 3 .W^. is the 

weight at Gauss point (r^s^). 


For 3-Gauss point quadrature in one direction 


s 1’ s 3 

= r 1 , r 3 = + 0.77459 

(4.30) 

s 2 = 

r 2 = 0.00000 

(4.31) 

*1 " 

W 3 = 0.55555 

(4.32) 

*2 - 

0.88888 

(4.33) 


4.4 EVALUATION OF BOUNDARY INTEGRALS: 

Referring to the boundary element shown in Fig 0 (4.4^ 

The x co-ordinate of any point on the line element 

is given by 

x = r 1 x 1 + r 3 x 3 (4.34) 

where, 

r 1 " 1 ”* 7 • ( 4 * 35 ) 

_ _ , x , 

r 2 - -JT! 


(4.36) 
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X 1 fXg the co-ordinates of extreme nodes 
^ ~ length of element. 

Then, 

m ' n ‘ 

, job n , m • n . ( e \ 

/ r. r 2 dx = . 

° ( m+n+1 ) \ 

Shape functions for the element are 

N 1 = r 1 ( 2r^ - 1) 

N 2 = 4^ r 2 

N 3 = r 2 ^ 2r 2 “ 1 ) 


(4.37) 


(4.38a) 

(4.38b) 

(4.38c) 


Using the results of eqns.(4.37 - 4.38c), the boundary 
integrals in the eqns. (4.15), (4.18b) and (4.18c) can be 

found as: 



(4.39) 

(4.40) 

(4.41) 


4.5 SOLUTION OF SIMULATIONS LINEAR DIFFERENTIAL EQUATIONS: 

To solve eqn . (4.14) for nodal temperature finite 
difference method is used for time domain,# 
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Taking part of forward and part of backward time 
derivative temperature at (n+l) th level and n th level are 
related by 


{T n i 


n+1 


3T n+1 

9 At l 


= n n r 


,n 


3T 

+ (1-9) At 


(4.42) 


where , 

A t = time step 

0 = 1, for fully backward or implicit scheme 
= 0 # for fully forward or explicit scheme 
= for semiimplicit scheme. 


Multiplying eqn. (4.42) by M n+1 and using implicit scheme, it 
can be written asi 


- 9T„ I n+1 _ +1 n 

[MJ n+1 JT n l n 1 - At [M] £ ^ ) | - [M] n tT n } n 

(4.43) 

From eqn. (4.14 ), 

[Mj n+1 t|^i n+1 = -[HliT a }| n+1 (4.44) 

From eqns . (4,43) and (4.44) 

([M] n+1 + At [ H) n+1 ) 5T n ! n+1 = [M] n+1 (T/ 

n+1 


+ At IF! 


(4.45) 
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Thus, knowing nodal temperatures at n th level it can 
be updated at n+1 level. Using eqn 0 (4,45), 

4.6 AVERAGE TEMPERATURE ALONG TOOL-CHIP INTERFACE: 

Envoking the parabolic variations along a 3— noded 
line element (Fig, (4, 4)), average temperature along tool- 
chip interface can be written as 

T ". “ 4 e ii nt (C Tdl (4 ' 46) 


whe re. 


H = chip contact at length 
]=int = number of line elements along interface 

L^ e ^ = length of element (e) . 


as 


Using the results of eqn, (4.40), we can write eqn,(4.46) 



E int 

£ 

e=1 


L (e) 

* 6 


( 






(4.47) 


where T 1 , T 2 , T 3 are the temperatures at local nodes 1 , 2 
and 3 respectively for the element (e)y> 



CHAPTER V 


RESULTS AND DISCUS5I ON 


A general purpose finite element programme has been 
developed to determine the temperature distribution in 
tool-chip-work system and average tool-chip interface temper- 
ature under steady state condition (i.e of longitudinal turning) 
as well as transient situation (i 0 e., accelerated cutting) 
for a wide range of cutting conditions. The results obtained 
have been compared with experimental [l2j and analytical [9] 
results to test the validity of present model „ The experi- 
mental cutting data are borrowed from the work of other 
researchers (the details are given in Appendix I). The vari- 
ation in temperature within the tool chip- workpiece system is 
illustrated by drawing isotherms. In order to investigate 
the effects of shear strain acceleration and its governing 
parameters on mean tool-chip interface temperature and temper- 
ature distribution in tool-chip-work system, accelerated 
test results have been compared with that of longitudinal 
turning. 

5.1.1 COMPARISON OF PRESENT MODEL WITH OTHER ANALYTICAL 

MODEL: 

Table 5.1 shows the comparison of the maximum temper- 
ature obtained using present analytical model and the 
results obtained by Murarka et al. [9 ] . It is clear from 



TABLE 5.1 : Comparison of Maximum Temperature Obtained, Using Present Model and Other 

Analytical Model during Longitudinal Turning 
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the table that the difference in temperature obtained by 
these two model ranges from 35°C to 58°C (about 8-1 C%), 
but the trend seems to be the same. Further temperature 
obtained using present work are always less than those 
obtained by Murarka [9 ] . This discrepancy could be 
attributed to the difference in the formulation of the 
problem. Murarka determined the velocity and strain rates 
experimentally using grid technique. The calculated 
tool-chip plastic contact length using experimental data. 

The present model does not require any of these experimental 
data. They have been calculated as described in Chapter 3. 

The value of maximum temperature is also governed by the 
fineness of the mesh in the critical zone (around the maximum 
temperature). "TV->e vesulir is also ccmp^ed ©la-Vajmed 

^ O.S. *slntf\OT> JTl ^ 5)0 - P ' T be ts ‘ 

^OtbdI eno > 

5.1.2 COMPARISON WITH EXPERIMENTAL RESULTS: 

The variation of the average tool-chip interface 
temperature (both analytical and experimental [123) with 
cutting speed for longitudinal turning and facing test are 
shown in Figs. (5.1a) and (5.1b) respectively. In longitudinal 
turning, the experimental and theoretical results are quite 
close to each other. The difference in temperature values 
ranges approximately from 5 to 15%. In facing, the temperature 
obtained theoretically is lower than the experimental values 
by an amount equal to about 15 to 25%. From these figures, 
it is clear that both the curves, experimental and theoretical, 
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show the same trend, i e e # an increase in the average tool- 
chip interface temperature with increasing cutting speed. 
However, a little discrepancy between the two exists. The 
difference between the experimental and theoretical values 
of temperature may be attributed to the following: 

(i) Since actual metal cutting process is very complex, 
several idealizations and approximations have been 
made. For example, tool is perfectly sharp, elastic 
contact length is neglected, two-dimensional model 
accurately explains the deformations and heat transfer 
etc , 

(ii) Experimental temperature results and cutting component 
of force (for cutting conditions shown in the Figs, (5,1a) 
and (5,1b) and depth of cut as 1,0 mm) have been 
borrowed from Reference [ 1 2 1, Feed force and shear 
angle data are not given in the above reference. They 
have been determined based on assumptions that the 
coeff icientof friction (H ) and shear angle (0) are same 
for depth of cut as 0,6 mm and 1,0 mm. The necessary 
data for depth of cut as 0 o 6 mm are availed from 
Reference [ 1 7 1 , other conditions being identical. These 
assumptions are solely responsible for discrepancy in 
the theoretical and experimental results, 

(lii) The actual machining process is dynamic, so statically 
calibrated thermocouple may not reflect the true value 



experimental experimental 

1200 ~~ theoretical ^ 1200 — theoretical 

f =0.26 mm/rev o f = 0. 26 mm/re 

N =100 rpm S N =100 rpm 


\ 

\ 

• \* 

\ 
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3 0 'clijuei aoDjje;ui diqo ioo; oOdjoav 


Cutting speed, m/s Cutting speed, m/s 

Fig. 5.1 Comparison of analytical and experimental average temperature 
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of temperature. Also, formations of BUE and oxide 
layers at the interface may yield misleading results. 
Thus, the use of thermocouple method for measuring 
the tool-chip interface temperature is not fully 
reliable 0 

5.2 TEMPERATURE DISTRIBUTION: 

The variation in temperature within the tool, chip 
and workpiece obtained using FE model is shown by isotherms 
for taper turning, facing and longitudinal turning in 
Figs„ 5.2a, 5.2b and 5.3 respectively. These results were 
obtained for cutting conditions shown in the figures. 

The isotherms in Figs. (5.2, 5.3) show that high 
temperature regions are the regions in the vicinity of the 
tool-chip interface. The lower temperature isotherms are 
located in the regions away from the interface. It may be 
explained as follows. High temperature regions are located 
near tool-chip interface because of the fact that the sum of 
frictional teat generated at tool-chip interface and heat 
generated due to chip deformation in secondary deformation 
zone is much higher than that of heat generated due to plastic 
deformation in PSDZ. The other reason for low temperature 
isotherms in and around PSDZ is that the large amount of 
heat generated in PSDZ due to plastic deformation is carried 
away by chip, and remaining is conducted in tool and work. 




Fig. 5.2(a) Temperature distribution in taper turning 
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Figures 5.4(a) and 5.4(b) show the variation of 
temperature during longitudinal turning, taper turning and 
facing along rake and flank face. It is evident from the 
figures that in the beginning, the temperature profile 
increases steeply with the increase in distance from the 
cutting edge along rate face. It attains a maximum value 
at a distance of about 0,6-0. 7 times the tool-chip contact 
length, thereafter it starts decreasing. The steep increase 
near the cutting edge is due to high heat generation rate 
in the secondary zone. But, beyond the maximum temperature 
point, the heat generation rate is smaller as compared to 
rate of heat loss due to conduction. Therefore, the temperature 
decreases with the increase in distance away from the maxima. 

The secondary deformation zone ends at a point where the tool 
and chip contact ends. Beyond this point, there is no heat 
generation, but there is heat loss due to conduction, and 
convection to the atmosphere. Therefore, the temperature goes 
on decreasing with the increase in distance. 

Almost similar pattern of variation in temperature 
along flank face is also observed. However, the magnitude 
of highest temperature and slope of temperature profile 
slightly differ in the two cases. The increase in temperature 
may be attributed to the fact that heat generation in the 
secondary zone also affects flank face temperature near the 
cutting edge. Also, there is some amount of heat generation 
due to friction between flank face and machined surface, 
however this has not been taken into account in the present model. 
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Then a fall of temperature is expected as one moves away 
from the cutting edge as a result of the heat loss due to 
conduction and convection. 

The similar trends (Fig, 5.4) of variation in temper- 
ature along rake and flank face have been reported by 
others C 9 , 111. In both the figures 5.4(a) and 5.4(b), it 
is also observed that the temperature attained at a distance 
from the cutting edge is highest in case of longitudinal 
turning lowest in case of facing and in between the two 
(longitudinal and facing) in case of taper turning. This 
conclusion is in line with the experimental observations 
reported in reference [12] and [20] about the average tool- 

chip interface temperature . “Tc*v, Was Wsen ealcu) eeFsd 

ei.eYYie.rrts amJ 2.2_ eJ evYve*vhs voe-sH, TVie. db-ffeaen co. ? -H/\£ 

ua 'tvOD C3!4.<ts o^eruz ~Y%* ven^- swvaM . 

5.3 PARAMETRIC STUDY IN LONGITUDINAL TURNING AND 
ACCELERATED CUTTING: 

To study the effects of strain acceleration governing 
parameters on mean tool-chip interface temperature during 
longitudinal turning and accelerated cutting, parametric 

studies has been made. 

Variation of average tool-chip interface temperature 
with feed rate and cutting speed, for longitudinal turning 
as well as accelerated cutting (viz. taper-turning and facing) 
are shorn in Figs. 5.5a and 5.5b respectively. It is clear 
from the figures that the average temperature in all the 
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Three cases, increases with the increase in feed rate and 
cutting speed. However, the magnitude of temperature is 
seen to be the highest in case of longitudinal turning 
and lowest in case of facing at any machining condi tion 0 
The similar trends have been observed by the investigators 
[12] and [20] as shown in Figs. 5.6a and 5.6b respectively. ■* 

Figure 5.5c shows the variation in average tool-chip 
interface temperature with taper angle in case of taper 
turning. The longitudinal turning and facing temperature 
results for the same cutting conditions (except the taper 
angle) are also shown in the Fig. 5.5c. It is clear from 
the figure that the temperature decreases with the increase 
in taper angle. 

The variation of average tool-chip interface temper- 
ature with spindle speed is shown in Fig. 5,5d. It is 
observed that the temperature increases with an increase in 
spindle speed in case of accelerated cutting* 

-X Tine miuTfi tsP vcxHjafU avr crP \~ri 

-t-ooft- oVd p ctce. 'is pmcnoT) "\&€ pccx&kblTc 

.spW C in -f^Tn i V* J -ftrr 

\ s. "not 'Sirn'i'icur. . *TV><® poss'jJole ^-a^sxAorx -fb®- oO^c-repa/rc^ 
ccnM Joe jtvre -CWrel -Prc^o V^-£>rJ . 
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Fig. 5.5 Comparison of tool -chip interface temperature during longi- 
tudinal turning, taper turning and facing. 
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CHAPTER 6 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 

6.1 CONCLUSIONS: 

From the results obtained in the present work, the 

following conclusions can be drawn: 

(1) Present model has been tested against analytical 
results for longitudinal turning and experimental 
results for longitudinal turning and facing test. 

It has been found that the correlation between the 
two is quite satisfactory. 

(2) It is observed that the pattern of variation of 
temperature, along rake and flank face of a tool 
during longitudinal turning, taper turning and 
facing is similar. Maximum temperature on the 
rake and flank face in all the three cases occurs 
at some distance away from the cutting edge. 

(3) It is also seen that mean tool-chip interface 
temperature is higher in case of longitudinal turning 
and lower in case of facing as compared to the temper- 
ature in taper turning for the same cuttin g conditions. 
Average tool-chip interface temperature in all the 
three cases are found to increase with the increase in 
cutting speed and feed rate. Increase in spindle 
speed results in an increase in average temperature 

in case of accelerated cutting. However, effects 
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of an increase in taper angle is to decrease the 
temperature in case of taper turning* 

(4) From the above, it can be concluded that the average 
tool-chip interface temperature, as well as the local 
value of temperature at any location within the tool- 
chip-work system are not same during longitudinal 
turning, taper turning and facing operation under the 
same set of machining conditions,* 

6.2 SUGGESTIONS FOR FUTURE WORK: 

(l) In the present analysis, the variation of physical 
properties (such as specific heat, density, thermal 
conductivity etc.) of tool and work materials with 
increase in temperature is not taken into account. 

This can be taken into to have better correlation with 
experimental results. 

(?) Experiments may be conducted for a varieties of tool- 
work material combinations in a wide range of cutting 
conditions. Using these experimental results, the 
accuracy of the present model may be established. 

(3) In the present analysis yield shear stress is assumed 
to be constant throughout in the PSDZ. In actual 
practice, it is not so. Hence a FE model can be 
developed which can take into account the variation 
in shear stress in PSDZ. This may yield more 
accurate results 0 
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(4) Effects of temperature on microhardness of chips, 
coefficient of friction along rake face, strain 
rate and strain acceleration should be studied 
in depth. 
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APPENDIX I 


EXPERIMENTAL DATA 


To test the validity of the analytical model developed 
in this text, experimental results have been borrowed from 
the work of Gupta [17] , Jain [12] and Murarka [ 9 ]. 

Using the experimental results, response surface models 
for shear angle, cutting force, feed force and shear flow 
stress were developed by Gupta [17], These models have been 
used to conduct parametric study given in Chapter 5. 

Jain [12] has reported the results of cutting force 
and average tool-chip interface temperature during longitudinal 
turning and facing. The tool and work materials and the machin- 
ing conditions in case of Gupta [ 1 7 ] and Jain [12] were almost 
same except the depth of cut. 

Temperature results in reference [12] is in terms of 
mV of e.m.f., which has been converted into degree centigrade 
using calibration chart (Fig, A1), taken from the reference [ 4 ], 
Experiments were conducted with d.c,c, as 1 mm. Due to 
unavailability of shear angle and feed force data [12], 
assumptions are made that the coefficient of friction (P ) and 
shear angle (0) remain the same for d,o,c, as 0,6 mm C 1 7 3 and 
d.o.c, as 1,0 mm [l 2 ] . Required data for 0,6 mm d.o e c. case 
are taken from ref,[l7]. 
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DESIGN OF EXPERIMENTS 

The design of experiment is the procedure of selecting 
the number of trials and conditions for running them. Number 
of trials should be sufficient for solving the problem that 
has been set with required precision. The Important features 
are, minimum total number of trials, simultaneous variations 
in all the parameters and the selection of clear cut strategy, 
permitting the investigator to make substantial decisions 
after each series of trials. 

To determine the linear, higher order and interaction 
effects of different parameters on shear angle, shearflow 
stress, cutting and feed force, experiments were conducted for 
the trials shown in Tables A1 , A2 and A3. The responses were 
measured or calculated and are given in Tables (A4, A5, A6). 

The second order response surface equation was fitted in the 
form 


B. 


n 

+ S 
i=1 


B i X i + 


n 

S 

i=1 


B it x t + i z <j b ±5 


x i Y j 


(A.1 .1) 


where the values of B’s are given in Table (A7) for temperature 
response. The labelled values of X^s are not used here (i.e. 
for temperature). 

* Temperature response included in the table is obtained using 
present finite element model. 
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LEVELS FOR DIFFERENT FACTORS 


TABLE A1 : Longitudinal Turning 


Levels 

Factors 

-2 

-1 

0 

1 

2 

Feed X 1 
( nun/ rev .) 

0.10 

0.15 

0.20 

0.25 

0.30 

Cutting X 2 

Speed 
(m/s ) 

0.34 

0.40 

0.46 

0.52 

0.58 

* Calculated values 

which 

were not 

available 

in the 



M/c Tool. 
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TABLE A2 : Taper Turning 


Levels 

Factors 

-2 

-1 

0 

1 

2 

Spindle Speed 
(R.P.M.) X, 

80 

160 

(240)*250 

320 

400 

Taper 

Angle x 

(Degrees) 

2° 

3° 

4° 

5° 

6° 

Feed Xg 

(ram/ rev) 

0.1 

0.15 

0.2 

0.25 

0.3 

Cutting X^ 

Speed (m/s) 

0.34 

0.4 

0.46 

0.52 

0.58 

Calculated values 

which 

were not 

available in 

the m/ c 

tool. 


81 


TABLE A3 : Facing 


i 

Factors 

Levels 

-2 

-1 

0 

1 

2 

Spindle 

(R.P.M.) 

Speed 

X 1 

80 

160 

(240)*250 

320 

400 

Feed 
( mm/ rev ) 

x 2 

0.1 

0.15 

0.2 

(0.25 )* 
0.26 

0.3 

Cutting 

Speed 

(m/s) 

X 3 

0.34 

0.4 

0.46 

0.52 

0.58 


* 

Calculated values which were not available in the m/c tool. 


PLAN OF EXPERIMENTS 
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TABLE A7 ! Values of Constants in Response Surface Model 
(Eqn. A. 1.1) for Average Tool-Chip Interface 
Tempe rature 


Constants 


Longitudinal Taper Facing 

turning turning 



B 24 

B 34 


-1529.3148 

1660.1075 

0.8214 

12.7846 

332.8345 

- 7820.2025 

9981 .9028 

-1 0709.4900 

-505 .41 96 

- 

- 0.0000 

0.0057 

14.4073 

10249.8070 

-503 .5757 

14535.1890 

6857 .3914 

- 

-1 .1 1 1 0 

-21 .2768 

-21 .2190 

-25.8536 

17.9569 

- 

1359.5124 

22819.4630 

-1 085.9270 

mm 

•21878.5380 

mm 


* Levelled values of factors X^’s in Eqn. (A .1.1) are not 
used. 
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